Strong color fields and heavy flavor production 



I. Bautista and C. Pajares* 
September 30, 2010 



Abstract 

The clustering of color sources provides a nat- 
ural framework for soft partonic interactions 
producing strong color fields. We study the 
consequences of these color fields in the pro- 
duction of heavy flavor and the behavior of 
the nuclear modification factor. 



1 Introduction 

Heavy flavor production in heavy ions colli- 
sions is an ideal probe to study the early time 
dynamics of these nuclear collisions. Several 
theoretical studies predict |lj a substantial 
enhancement of open charm production asso- 
ciated to deconfined parton matter relative to 
the case of a purely hadronic scenario with- 
out quark-gluon plasma formation. Recent 
studies point out that the dynamics of heavy 
quarks is dominated by partonic interactions 
in a strongly coupled plasma modeled neither 
by hadronic interactions nor by color screen- 
ing alone [2J. Therefore, these quarks are very 
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relevant in the study of the initial state of the 
collision. Owing to difficulties in reconstruct 
the D-mesons decay vertex, RHIC experi- 
ments have measured open charmed quarks 
indirectly, via the semileptonic decay to non- 
photonic electrons or muons [3] [I]. In the 
standard picture charm quarks are produced 
by initial gluon fusion and their produc- 
tion rates are expected to be well described 
by perturbative quantum chromodynamics 
(pQCD) at fixed order plus next-to-leading 
logarithms (FONLL) [5]. The suppression 
of single, nonphotonic electrons or muons at 
RHIC is usually attributed to heavy-quark 
energy loss. As a charmed quark of energy E 
cannot radiate gluons forming an angle below 
arcsin(m/i?) (dead cone effect), it is expected 
that heavy quarks lose less energy than light 
quarks [6], but the suppression experimen- 
tally observed is similar. In fact, many cal- 
culations based on energy loss via hard scat- 
tering [7] or via multiple soft collisions [8] ob- 
tained less suppression than the experimental 
data when the beauty contribution is taken 
into account. Similar results are obtained in 
evaluations based on medium interactions or 
collisional dissociation [9]. However, it has 
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been argued [TO] [EE] that under the assump- 
tion of an enhancement of the heavy-quark 
baryon-to-meson ratio, analogous to the case 
of the proton-to-pion and the A-to-kaon ra- 
tios measured in Au-Au collisions at RHIC, 
it is possible to achieve a larger suppres- 
sion of the nuclear modification factor. This 
is possible, because the heavy-quark mesons 
have a larger branching ratio to decay inclu- 
sively into electrons as compared to heavy- 
quark baryons and therefore, when the for- 
mer are less copiously produced in a heavy- 
ion environment, the nuclear modified fac- 
tor decreases, even in the absence of heavy- 
quark energy loss. Indeed the single nonpho- 
tonic nuclear modified factor, R e AA can be ex- 
pressed as [12] R e AA = R AA AF where R AA A is 

the nuclear modified factor for D and A c , i.e. 



pD+A 
n AA 



N% A + N\ A 



with N D and A^ A the produced D and A in 
AA or pp collisions and the N co u is the num- 
ber of collisions at a given centrality. The 
factor F is given by the expression 

= (l + a)(l + xCa) 
" (1 + Ca)(l + xa) { ' 

where, a and Ca are the charmed baryon-to- 
meson ratios in proton-proton and AA col- 
lisions respectively. Therefore C represents 
the enhancement factor for the ratio of charm 
baryons to mesons in AA as compared to pp 
collisions, x is the ratio between the branch- 
ing ratios for the inclusive decay of A and D 
into electrons: 



In [12] x has been estimated to be 0.14. As 
long as C is larger than 1 the factor F be- 



pD+A 



The 



comes less than 1 and R e AA < 
main question to solve is whether the ex- 
pected charmed baryon-to-meson expected 
enhancement is large enough to explain the 
difference with the experimental data. 

In a high-energy heavy-ion collision, strong 
color fields are expected to be produced be- 
tween the partons of the projectile and tar- 
get [13] [33] [IS] . These color fields are similar 
to those that appear in the glasma [16] pro- 
duced in the color glass condensate (CGC). 
In a string heavy-quark pairs are produced 
via the Schwinger mechanism with a rate 

Tqq = exp -rr- where k, is the effective 

string tension, proportional to the strength of 
the field (for a single string k\ ~ 1 GeV/fm). 
Longitudinal string models predict for heavy 
flavor a very suppressed production rate, 
since 



QQ 



qq 



exp i 



-7T 



(m. 



Q 



ml)) « 1 (4) 



(3) 



for Q = c and q = u,d. The color in these 
strings is confined to a small area in the trans- 
verse space, ttTq, with r — 0.25fm. In a 
central heavy-ion collision many strings are 
formed between the partons of the projectile 
and target in a limited collision start- 
ing to overlap each other, forming clusters. 
The field strength of the cluster is propor- 
tional to the square root of the number of 
strings. So, for a cluster of nine strings, the 
string tension increases more than eight or- 
ders of magnitude becoming comparable to 
the initial FONLLpQCD. The effect of strong 
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color fields for open charm has been investi- 
gated before [T3] showing that a three- fold 
increase of the effective string tension results 
in a sizable enhancement of the total charm 
cross section and the nuclear modified factor 
shows a suppression at moderate px consis- 
tent with the RHIC data. 

In this paper, we study the effects of strong 
color fields in the framework of percolation of 
strings [T7] . In this framework, a strong color 
field is obtained inside the clusters formed by 
the overlapping of individual strings. The 
clusters behave like individual strings with 
a higher string tension owing to the higher 
color field, and their energy momentum is 
the sum of the energy-momenta of individ- 
ual strings. The color field of the string is 
stretched between a quark and an antiquark 
or between a diquark and a quark located at 
the extremes of the string. In the case of a 
cluster, instead of quarks or antiquarks we 
have complexes Q and Q formed from the 
different quarks and their antiquarks or di- 
quarks and quarks at the extremes of the in- 
dividual strings [H][T9]. The clusters behave 
like strings with complexes QQ, located at 
the end, decaying into new pairs QQ, QQ, 
until they come to objects with mass com- 
parable to hadron masses, which are identi- 
fied with the known hadrons by combining 
the produced quarks or antiquarks with the 
appropriate statistical weights. In this way, 
the production of baryons and antibaryons is 
enhanced with the number of strings in the 
cluster. The cluster not only has a stronger 
color field than the individual string giving 
rise to a mass-enhancement effect but also 
enhances the production of baryons relative 



to mesons owing to the increasing probabil- 
ity of getting three quarks or three antiquarks 
from the complex QQ [18] . This second effect 
is similar to what happens in coalescence or 
recombination models [20] [21] . 

The percolation of strings incorporates to 
some extent the recombination of flavors in a 
dynamical way. Indeed a dynamical quark re- 
combination model has shown a sizable sup- 
pression factor for the nonphotonic electrons 
nuclear modification factor [12] . 

We evaluate the nuclear modification fac- 
tor for D Q , A c and B at RHIC energies, com- 
puting also the baryon-to-meson ratio in AA 
and pp collisions. We observe also in pp a 
moderate enhancement of the ratio func- 
tion of the transverse momentum which has 
consequences concerning the value of F and 
therefore the rate of the nonphotonic electron 
suppression. The plan of the paper is as fol- 
lows: In the next section, we introduce briefly 
the percolation of the strings, and then we 
present our results and conclusions. 

2 The string percolation 
model 

In the string percolation model [17j[22j[23j 
[21] [25], multiparticle production is de- 
scribed in terms of color strings stretched 
between the partons of the projectile and 
the target. With increasing energy and/or 
atomic number of the colliding particles, the 
number of strings N s , grows and they start 
to overlap forming clusters, very much like 
disks in two-dimensional percolation theory. 
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At a certain critical density, a macroscopical 
cluster appears, which marks the percolation 
phase transition. This density corresponds 
to the value 1] C = 1.2 — 1.5 (depending on the 
profile function of the colliding nuclei) where 
rj = N s Si/Sa and Sa stands for the over- 
lapping area of the colliding objects. A clus- 
ter of n strings behaves like a single string 
with the energy momentum corresponding to 
the sum of individual ones and with a higher 
color field corresponding to the vectorial sum 
in color space of the color fields of the individ- 
ual strings. In this way, the mean multiplicity 
< fi n > and the mean transverse momentum 
squared < p\ n > of the particles produced by 
a cluster are given by 



<l*i> 



(5) 



and < >= J ^ < v\\ >> where < fi\ > 
and < p T i > are the corresponding quantities 
in a single string. 

In the limit of high density of strings, equa- 
tions (5) transforms into [24] 



< fj, >= N s F(r)) < /ii > 



(6) 



<Pri> 



with F{rj) 



l-e—n 



F(v) v,/ y v 

For a specific kind of particle i, we will use: 
< fix >i, < p 2 T1 > u < fx n >i, and < p\ n >, 
for the corresponding quantities. To com- 
pute the multiplicities, we must know N s and 
fii (so for a fixed centrality, knowing N s we 
deduce the density rj). Up to RHIC ener- 
gies, in the central rapidity region N s is ap- 
proximately twice the number of collisions, 
Neon. However N s is larger than 2N co u at 



RHIC and LHC energies, in the same way 
as in nucleon-nucleon collisions. According 
to color exchange models, such as dual par- 
ton model or the quark gluon string model 
[26j[27], the number of produced strings N s 
is larger than two, starting at the RHIC 
energies. Indeed, at high enough energy 
the strings are stretched not only between 
the diquarks(quarks) and quarks(diquarks) 
of the projectile and target respectively, but 
also between quarks(antiquarks) and anti- 
quarks(quarks) of the sea. As the energy in- 
creases, more q — q or q — q are formed and 
N s becomes larger than two. For the same 
reason in AA collisions, N s at high energy is 
larger than 2N co u. In this work we take the 
values of N s from a Monte- Carlo based on the 
quark- gluon string model [2%] . 

Concerning the transverse momentum dis- 
tribution, one needs the distribution g(x,pr) 
for each cluster, and the mean square trans- 
verse momentum distribution of the clusters 
W(x), where x is the inverse of the mean 
of the squared transverse momentum of each 
cluster which is related to the cluster size by 
equation (5). We take g(x,p^) = exp (— p\x) 
as it is used for fragmentation of the Lund 
string. For the weight function W(x) we 
take the gamma distribution. The general- 
ized gamma distributions are unique distri- 
butions stable under the cluster-size transfor- 
mations |22j[29j[30j; for simplicity we choose 
gamma distribution simplicity 



with 



W(x) = 1 exp (-kx) (7) 



r(*) 



7 = hj < x > 



(8) 
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and 



1 

k 



< x 2 > - < x > 2 



(9) 

< x > z v ' 

The function k measures the width of the 
distribution (7) and is the inverse of the nor- 
malized dispersion of the transverse momen- 
tum squared. The function k depends on the 
density of strings r\. 

The transverse momentum distribution 
f(pr,y) of particle i is 



f(pr,y) 



dN 
dp\dy 



dxW(x)g(pT, x) 



(10) 



dN /^l T?( \ 1 

k r VD (l+F(r,)pl/k< P ^> u ) k 



dy 



The formula (10) is valid for all types of 
collisions, energies and also all kind of fla- 
vors. Later we will extend (10) for baryons. 
The function k(r]) was determined by com- 
paring (10) to RHIC data. The function k 
decreases with rj up to values r] ~ 1 (periph- 
eral Au-Au collisions at RHIC energies) from 
there it increases slowly. This behavior was 
expected. In fact, at low density there is no 
overlapping of strings and there are isolated 
strings; therefore k — > oo. When the density 
and therefore the numerator of Eq. (9) in- 
creases then k decreases. The minimum of k 
will be reached where the fluctuations in the 
cluster-size reach its maximum. Above this 
point, increasing rj these fluctuations decrease 
and k increases. The agreement with data for 
Pt up to 5 GeV/c is very good [22j[23j. 

In percolation of strings the fragmenta- 
tion of a cluster of many strings is via the 
Schwinger mechanism, producing successive 
pairs QQ, where Q represents the complexes 



of quarks, diquarks and antiquarks at the 
extremes of the original string. It is clear 
that formula (10) only contains the effect of 
the stronger color field of the cluster, which 
enhances heavy particles production, irre- 
spective of their being mesons or baryons, 
but it does not contain the breaking via fla- 
vor complexes QQ and therefore cannot de- 
scribe baryons correctly. In previous pa- 
pers [18] [19] [28] Monte Carlo codes were pre- 
sented where this mechanism was built up, 
but with the approximation of fusion of only 
two strings [T7] [18] or using an effective color 
field |28j. To to keep a closed analytical 
formula, incorporating the antibaryon and 
baryon enhancement from the mechanism de- 
picted here, we observe that this enhance- 
ment is similar to using the formula (10) 
with a larger density, or equivalently with a 
larger N s . This means that for antibayons or 
baryons if we want to continue with formula 
(10) we must replace r\ by t/b, 

fj B = Nfr) (11) 

and instead of the first equation (6) we must 
use 

m = Nt +a F( VB )v 1B (12) 

where the parameter a is obtained from a fit 
to the experimental dependence of the px- 
integrated p spectra with centrality [31J. The 
obtained value is a = 0.09. In same sense we 
can say that the antibaryons (baryons) probe 
a higher density than mesons for a fixed en- 
ergy and type of collisions. 

The equations (10), (11) and (12) allow us 
to compute the antibaryon (baryon) spectra. 
The equations (11) and (12) replace the re- 
combination process described in this section 
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and it should be considered as an approxima- 
tion to keep the analytical formula (10). The 
formulae (10), (11) and (12) are valid for all 
kind of particles and not only for heavy fla- 
vor. We will show some results concerning 
light flavor. 

3 Results 

The equation (10) is limited to low and mod- 
erate pt not higher than 4 — 5 GeV/c. In 
fact, we consider a gaussian pt distribution 
for the particles produced from the fragmen- 
tation of a string, without any power-like tail. 
This excludes the high-p^ behavior, although 
our formula (10) allows for an interpolation 
from low to high p?. By continuity, the high 
Pt suppression observed at RHIC should give 
rise to a suppression at moderate px, say 4 — 5 
GeV/c, which is the limit where our equa- 
tions apply. 

To know the pt distributions given by for- 
mula (6) we need the values of < p T >id~< 
Pt >\ d and < p\ >ia c ^< Pt >? Ac , i. e. 
the mean pt of D and A c particles produced 
from one string. We use < pt >id— 1.5 
GeV/c and < p\ > lk = 1.9GeV/c. The 
difference between these two values is close 
to the difference between the masses of Dq 
and A c and also agrees with the difference 
between the values commonly used of pri- 
mordial transverse momentum of pions and 
protons, < pt >i-k= 0.2 — 0.3 GeV/c, < 
Vt >ip— 0.6 — 0.7 GeV/c. For B we use 
< PT >1B = 4.25 GeV/c. 

In formula (10) the normalization is es- 
tablished by the values of ^f- at pt = 0, 




Figure 1: Raa for Au+Au central collisions, 
bars is data taken from PHENIX [3T]. 

which are computed using the formulae (6) 
for D and (9) for A c . To do this, we use 

2 

the values u,id = exp(— F(rj) , m °i )(J>in and 

<PT>i£, 

/iiA c = exp{-F{r] Kc )m\J < p T >ia b W> 
with fii n = 0.8 [25]. We use these func- 
tions for \i\T> and /i\ c because for heavy par- 
ticles, m T is very different from p\. Con- 
cerning the function k(rj), we take the shape 
and values from the studies done in reference 
[25] for AA collisions. In the case of pp we 
take k(rj) = 3.97 at </s = 200 GeV and 
k(rj) = 4.07 at y/s = 5.5 TeV. We discuss 
later the sensibility of the obtained result for 
the ratio (A c / D°) to different k values. 

In fig. 1 we present our results for the 
nuclear modified factor Raa for Au-Au col- 
lisions at RHIC for D°, A c , B and R e AA using 
formulae (1) as a function of p T compared 
with the experimental data on nonphotonic 
electrons. 

The overall normalization is given by the 
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value of Raa at = which has to do with 

2 

the factor exp (F(r] pp ) - F(r) AA )) <p ?> lg ■ 
Since we know the number of strings pro- 
duced in pp and AA collisions, we know r] pp , 
r/AA, F(jiaa) and F(i] pp ) and the only free 
parameter is < p\ >m- From the data we 
obtain < pt >id~ 1-5 GeV/c. The experi- 
mental errors allow us a 15% freedom in the 
value of < pt >id, however a higher value 
than 1.5 GeV/c would not be realistic an a 
lower value will give rise to a higher normal- 
ization and therefore Raa for Pt > 4GeV/c 
will exceed the experimental data even more 
than with the used value. The value of Raa 
for D° at low p? agrees with the results in 
[14|[32j, for p T > 4GeV/c we obtain an R AA 
larger than the non-photonic leptonic data. 
In fig. 2, we present our results on Raa at 
yfs = 5.5 TeV for a D°, A c and B. We see 
that, as expected, as energy increases the low 
Pt Raa increases, although the suppression 
at intermediate pr is similar. 

In fig 3, we present the ratio A c / D° for Au- 
Au at y/s = 200GeV. We observe that the ra- 
tio increases up to a maximum of 1.45 around 
p T 4 — 5 GeV/c. A very similar enhance- 
ment has been obtained in the dynamical re- 
combination model 1 12 . 

For comparison we include also our results 
for p/ir at central Au-Au collisions together 
with experimental data [33]. In fig. 4 we 
show the ratio A c /D° for pp at y/s = 200 
GeV. We observe a very smooth enhance- 
ment. In fig 5 and fig. 6 we show our re- 
sults for the ratio A c / D° at y/s = 5.5TeV 
for Pb-Pb collisions and pp collisions respec- 
tively. We observe in both of them larger 




Figure 2: Raa for Pb+Pb central collisions 
at y/s = 5.5 TeV. 



1.4 — 
1.2 — 




Figure 3: Squares are used for ratio A c /D°, 
starts are used forp/71", and errorbars are used 
for data from PHENIX for Au-Au central col- 
lisions at v/i = 200 (GeV). 
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Figure 4: Ratio A c /D° for p+p collisions at 
y/s = 200 GeV. 



enhancement than at RHIC energies, partic- 
ularly in the nuclear case. In fig. 7 we plot 
the factor F at y/s = 200 GeV (blue) and at 
y/s = 5.5 TeV (red line). We observe that at 
RHIC energies the factor F is only slightly 
below one, and for — 4 — 5 GeV/c it is 
clearly over 0.5, which means that the A c /D° 
enhancement in AA is not able to explain 
all the difference between experimentally ob- 
served suppression of Raa for non-photonic 
electrons and the pQCD expectations. We 
have studied the effects due to the uncertain- 
ties in the k values for pp at this energy For 
reasonable alternative k values the enhance- 
ment of A c /Dq in pp with px is larger giv- 
ing rise to lower C factor in equation (3) and 
therefore the factor F is near to one, consis- 
tent with our main conclusion, namely that 
the A c / D° enhancement is not able to explain 
all the difference between the experimentally 
observed values and the perturbative expec- 
tations. 




Figure 5: Ratio A c /D° for Pb-Pb central col- 
lisions at y/s = 5.5 TeV. 




P T tGeV/c) 



Figure 6: Raa 
yfe = 5.5 TeV 



A c /D° for p+p collisions at 
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Figure 7: Factor F for central collisions at 
RHIC in full line and LHC energies in dashed 
lines. 

4 Conclusion 

The overlapping of the strings formed in the 
collision of heavy nuclei particles produces 
strong color fields which give rise to an en- 
hancement of heavy flavor. We have com- 
puted the nuclear modification factor of D°, 
A c and B° at RHIC and LHC energies for 
AA collisions. Referring to D°, we obtain a 
good agreement at low p T with the experi- 
mental data for the nuclear modification fac- 
tor of non-photonic electrons. For px values 
between 2 and 6 GeV/c our results obtained 
are over the experimental data as pQCD |18J . 

The ratio A c / D° as a function of pt for Au- 
Au collisions at y^i = 200 GeV is enhanced 
showing a maximum around 5 GeV/c. Such 
as enhancement is much larger at LHC ener- 
gies. However, the enhancement A c /D° can 
explain only half of factor 2 difference be- 
tween the experimental data and the pQCD 



expectations at RHIC energies. 

In p-p collisions the ratio A c /D° also rises 
SIS db function of px but very smoothly at 
RHIC energies. At LHC this increase is a fac- 
tor of 2 between px = and px = 6 GeV/c. 
The enhancements of A c /D° in A A and pp 
collisions are larger at LHC than at RHIC 
as it was expected due to the stronger color 
fields produced. 
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